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2.  Status  of  Effort 

The  goals  of  this  research  program  are  to  provide  a  detailed  understanding  of  the 
recombination  mechanism  of  atmospheric  pressure  air  plasmas  and  to  assess  the  reaction 
rates  proposed  in  the  literature.  During  previous  years  of  this  program,  it  was  determined 
that  overall  plasma  recombination  is  mostly  controled  by  the  rate  of  a  single  reaction,  the 
three-body  recombination  of  NO.  Experiments  with  a  recombining  air/argon  plasma  were 
then  designed  in  order  to  provide  an  assessment  of  this  key  rate.  Several  optical  diagnos¬ 
tic  techniques  were  developed  to  interrogate  important  plasma  parameters  including  tem¬ 
peratures  and  species  concentrations.  In  particular,  a  new  technique  of  electron  number 
density  measurements  based  on  measurements  of  the  nonequilibrium  population  of  the 
predissociative  C  state  of  NO  was  devised.  The  results  obtained  over  the  past  year  are 
presented  in  this  report.  It  was  found  that  the  three-body  NO  recombination  rate  pro¬ 
posed  by  Dunn  and  Kang  appears  to  be  too  slow  by  a  factor  -100,  and  that  the  rate  pro¬ 
posed  by  Gupta  et  al.  or  Park  is  also  too  slow,  by  about  a  factor  -10.  This  result  has  sig¬ 
nificant  implications  for  the  development  of  techniques  to  create  and  maintain  elevated 
electron  number  densities  in  air  plasmas. 


3.  Accomplishments/New  Findings 


3.1.  Background 

As  discussed  in  our  previous  progress  reports,  electron  recombination  in  atmos¬ 
pheric  pressure  air  plasmas  at  temperatures  between  4,500  and  7,500  K  occurs  primarily 
through  fast  and  equilibrated  two-body  dissociative  recombination  reactions,  mainly: 

NO'"  +  e  o  N O,  (1) 

and  therefore  the  extent  of  ionizational  nonequilibrium  under  these  conditions  is  gov¬ 
erned  by  the  rates  of  the  (slow)  neutral  recombination  reaction, 

N  +  O  +  M  NO  -I-  M.  (2) 

Thus,  the  rate  of  Reaction  (2)  determines  the  extent  of  ionizational  nonequilibrium  in  air 
plasmas.  Yet,  very  large  discrepancies  exist  in  the  rate  proposed  for  this  reaction  by  vari¬ 
ous  authors,  namely  Dunn  and  Kang,*  Gupta  et  al.,^  and  Park.^’^  The  rate  proposed  by 
Park  is  somewhat  greater  than  the  rate  of  Gupta  et  al.,  which  is  itself  greater  than  Dunn 
and  Kang’s  rate  by  a  factor  10. 

To  assess  the  predictions  of  three  air  kinetics  mechanisms,*-^  experiments  were 
conducted  with  recombining  atmospheric  pressure  air  plasmas  generated  with  the  50  kW 
inductively  coupled  radio  frequency  plasma  torch  described  in  our  previous  progress  re¬ 
ports.  Initially  in  equilibrium  at  a  temperature  of  7,160  K,  the  plasma  was  cooled  to 
-4,900  K  within  -400  )ts  by  flowing  at  velocities  approaching  0.5  km/s  through  water- 
cooled,  brass  test-sections  mounted  on  the  plasma  torch  nozzle  exit.  Through  detailed 
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measurements  of  temperatures,  electron  number  densities,  and  emission  spectra,  it  was 
possible  to  conclude  that  the  rate  of  Reaction  2  should  be  faster  than  the  rates  proposed 
by  Dunn  and  Kang  and  Gupta  et  al,  and  at  least  as  fast  as  the  rate  proposed  by  Park. 

To  provide  a  better  assessment  of  the  rate  of  Reaction  2  (on  which  the  pure  air 
studies  permitted  us  to  place  a  lower  limit),  and  to  extend  the  range  of  our  previous  work 
to  lower  temperatures,  experiments  were  conducted  with  a  -10%  air  -  90  %  argon  plasma 
(15.6  slpm  air,  162  slpm  argon,  and  2.3  slpm  H2)  cooled  from  7,900  K  to  approximately 
2,500  K  over  a  time  of  -1.3  ms  and  a  distance  of  65  cm.  The  small  quantity  of  H2  was 
premixed  into  the  gas  for  purposes  of  electron  number  density  measurements.  Dilution  of 
-10%  air  in  argon  causes  the  recombination  of  nitrogen  atoms  to  be  about  one  order  of 
magnitude  slower  than  in  the  pure  air  plasma  case  as  the  third  body  efficiency  of  argon  in 
Reaction  2  is  approximately  20  times  smaller  than  the  third  body  efficiencies  of  N  and  O 
atoms.  By  slowing  the  overall  N  atom  recombination,  it  should  be  possible  to  observe 
significant  chemical  and  ionizational  nonequilibrium,  and  thus  to  place  an  upper  limit  on 
the  rate  of  Reaction  2.  Even  at  this  relatively  high  dilution  level,  the  electron  recombina¬ 
tion  path  remains  the  same  as  in  pure  air  because  electrons  still  recombine  preferentially 
via  the  dissociative  recombination  reaction  NO^  +  e  ^  N  +  O,  even  though  NO^  is  no 
longer  the  dominant  ion  here. 


3.2.  Air/Argon  recombination  experiments 

In  this  section,  the  results  of  our  experimental  studies  with  an  air/argon  plasma  are 
presented  and  discussed.  For  these  measurements,  optical  emission  was  collected  at  the 
exit  of  test  sections  ranging  from  0  cm  (nozzle  exit)  to  65  cm  in  length.  As  the  tempera¬ 
ture  varies  from  -8000  K  at  the  nozzle  exit  to  -2500  K  at  the  65  cm  test  section  exit,  sev¬ 
eral  different  techniques  were  developed  to  obtain  temperatures  and  electron  number  den¬ 
sities  at  these  locations. 

0.  10.  and  15  cm. 

At  these  locations,  the  plasma  temperature  was  sufficiently  high  to  obtain  tem¬ 
perature  profiles  from  the  Abel-inverted  line  intensities  of  oxygen  (777.3  nm),  argon 
(763.5  nm)  and  hydrogen  (Ha).  At  each  location,  the  temperatures  determined  from  the 
three  atomic  lines  were  found  to  agree  within  150  K.  The  temperature  profiles  measured 
at  0,  10  and  15  cm  are  shown  in  Fig.  1. 
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Figure  1.  Temperature  profiles  of  the  air/argon  plasma  at  0,  10,  and  15  cm  downstream  of  the  nozzle  exit. 


The  centerline  temperature  at  these  locations  is  also  sufficiently  high  that  a  spec¬ 
trally  resolved  measurement  of  the  Hp  line  shape,  and  thereby  of  the  electron  number  den¬ 
sity, ^  can  be  obtained.  The  measured  electron  number  densities  match  the  LTE  number 
densities  (Fig.  2)  based  on  the  measured  centerline  temperatures  at  0  and  10  cm.  At  15 
cm,  an  electron  overpopulation  of  ~3  is  measured. 
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Figure  2.  Equilibrium  and  measured  electron  number  densities  (air/argon  plasma). 


Figure  3  shows  comparisons  between  the  equilibrium  line-of-sight  spectral  simu¬ 
lations  based  on  measured  temperature  profiles  and  the  measured  line-of-sight  emission 
spectra.  The  high  centerline  temperatures  at  these  locations  necessitate  consideration  of 
the  small  but  noticeable  contribution  of  the  electron  recombination  continuum  to  the 
measured  spectra.  Due  to  the  uncertainties  of  accurately  calculating  this  continuum  over 
the  measured  temperature  and  wavelength  regimes,  the  continuum  is  modeled  with  a  con¬ 
stant  value  added  to  the  NEQAIR2  simulations  such  that  the  computations  and  measure¬ 
ments  match  at  360  nm.  These  small  additive  constants  are  10,  0.7,  and  0.1  p,W/cm^/sr  at 
the  0,  10,  and  15  cm  locations,  respectively.  As  discussed  in  Appendix  A,  these  values 
are  consistent  with  estimates  of  the  argon  recombination  continuum. 
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Wavelength  (nm) 

Figure  3.  Measured  and  computed  LTE  emission  spectra  of  the  air/argon  plasma  at  0,  10  and  15  cm. 

As  can  be  seen  in  Fig.  3,  the  line-of-sight  emission  spectra  at  0,  10,  and  15  cm  are 
reproduced  within  the  20%  experimental  uncertainty  by  equilibrium  spectral  simulations 
at  the  measured  temperatures  at  almost  all  wavelengths.  At  0  cm,  the  disparity  between 
the  two  curves  near  200  nm  is  attributed  to  the  inaccuracy  in  the  assumption  of  a  constant 
electron  recombination  continuum.  Both  the  argon  continuum  model  and  the  empirically 
determined  additive  constants  show  a  strong  dependence  on  electron  number  density  and 
a  weaker  dependence  on  temperature.  As  the  electron  number  density  and  temperature 
decrease  with  distance  from  the  nozzle  exit,  the  magnitude  of  the  continuum  becomes 
negligible  in  comparison  with  molecular  radiation. 
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At  all  three  locations,  the  absolute  intensity  of  the  OH  A  X  ^Hi  band 

(280<X,<  315  nm)  is  not  well  reproduced  by  the  spectral  simulations.  This  discrepancy 
can  be  explained  by  considering  the  LTE  number  density  of  the  OH  A  state  (Fig.  4)  which 
has  a  maximum  at  approximately  3700  K,  well  below  all  temperatures  at  a  radius  of  0.35 
cm.  The  discrepancy  therefore  arises  as  the  measured  OH  emission  at  the  0,  10,  and  15 
cm  locations  originates  in  the  outer  edges  (r  >  0.4  cm)  of  the  plasma  where  the  Abel- 
inverted  atomic  line  intensities  are  too  weak  to  provide  reliable  temperature  information 
for  the  spectral  simulations. 

In  summary,  at  0  and  10  cm  all  measurements  are  consistent  with  LTE.  At  15  cm, 
the  spectral  emission  is  also  consistent  with  LTE  despite  the  small  measured  electron 
overpopulation. 


Temperature  (K) 

Figure  4.  LTE  OH  A  state  concentration  for  the  air/argon  plasma. 

40.  50.  65  cm 

At  these  locations,  the  excited  state  populations  of  atomic  species  are  so  small  that 
atomic  lines  cannot  be  seen  in  emission.  However,  the  temperatures  can  still  be  accu¬ 
rately  deduced  by  examining  the  temperature  sensitive  shape  of  the  measured  OH 
A^Z'^^X^Hi  band.  In  particular,  the  relative  strengths  of  the  two  peaks  at  306.8  and 
309.2  nm  strongly  depend  on  rotational  temperature  as  both  peaks  belong  to  the  (0,0) 
band.  The  decrease  in  intensity  beyond  311  nm  is  sensitive  to  both  vibrational  and  rota¬ 
tional  temperatures.  As  the  centerline  temperature  for  these  locations  was  found  to  be 
less  than  3700  K,  the  modeling  of  the  measured  line-of-sight  OH  A^Z""^  X^H  emission 
spectra  does  not  encounter  the  difficulties  discussed  previously. 

The  centerline  temperature  is  determined  by  comparing  numerical  OH  spectra^  at 
various  temperatures  normalized  to  the  309.2  nm  peak  to  the  measured  line  of  sight  spec¬ 
tra.  This  shape-matching  technique  yields  centerline  temperatures  with  an  accuracy  of 
-200  K.  The  remainder  of  the  temperature  profile  is  obtained  from  the  Abel-inverted  in¬ 
tensity  of  the  306.8  peak  relative  to  the  304  nm  baseline.  The  measured  line-of-sight  OH 
spectrum  at  40  cm  is  compared  with  the  synthetic  spectrum  in  Fig.  5.  The  agreement 
between  the  absolute  calculated  and  measured  intensities  is  typical  of  all  locations. 
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Wavelength  (nm) 

Figure  5.  Comparison  of  measured  and  calculated  line-of-sight  OH  A— >X  band  at  40  cm  (air/argon 
plasma). 


At  50  and  65  cm  this  procedure  must  be  slightly  modified  as  the  plasma  was 
found  to  be  in  chemical  and  ionizational  nonequilibrium.  The  shape  of  the  measured 
spectra  still  determines  the  centerline  temperature.  However,  an  overpopulation  factor  for 
the  OH  A  state  (poH  a)  must  be  determined  such  that  the  numerical  and  measured  line-of 
site  spectra  agree  absolutely.  This  method,  of  course,  assumes  that  the  overpopulation 
factor  is  uniform  across  the  plasma.  The  OH  temperature  profiles  so  determined  are  dis¬ 
played  in  Fig.  6. 


Radial  Position  (cm) 

Figure  6,  Temperature  profiles  of  the  air/argon  plasma  at  40,  50,  and  65  cm  downstream  of  nozzle  exit. 
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Intensity  (|LiW/cm  sr)  Intensity  (|iW/cm  sr)  Intensity  (|LiW/cm  sr) 


Features  emanating  from  the  NO  A,  C,  and  D  states  and  from  the  N2  C  state  are 
also  visible  at  these  locations.  Their  shapes  are  not  as  sensitive  to  rotational  or  vibra¬ 
tional  temperature  as  the  OH  spectrum  but  can  be  used  to  confirm  the  deduced  tempera¬ 
tures  within  500  K.  Electronic  overpopulation  factors  are  determined  in  the  same  manner 
as  for  the  A  state  of  OH. 


200  210  220  230  240  250  260  270  280  320  330  340  350  360  370  380  390  400 


200  210  220  230  240  250  260  270  280  320  330  340  350  360  370  380  390  400 


200  210  220  230  240  250  260  270  280  320  330  340  350  360  370  380  390  400 


Wavelength  (nm)  Wavelength  (nm) 

Figure  15.  Measured  and  computed  emission  spectra  of  the  air/argon  plasma  at  40,  50,  and  65  cm.  Included  in  the 
spectral  simulations  are  the  NO  p  (B— >X),  y  (A^X),  8  (C— >X),  e  (D^X),  O2  Schumann-Runge  (B— >X),  N2  2* 
(C->B),  Nj  r(B^X),  OH  (A->X),  and  NH  (A->X)  transitions.  Unless  specified  otherwise,  the  nonequilibrium 
factors  for  all  species  are  set  equal  to  1  in  the  simulations. 
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The  experimental  (gray  line)  and  numerical  line-of-sight  spectra  at  the  40,  50  and 
65  cm  locations  are  shown  in  Fig.  7.  The  region  between  280  and  320  nm  is  excluded  as 
the  much  more  intense  OH  A-X  band  (13,  8,  and  3  |iW/cm^/sr  respectively)  dwarves  the 
other  spectral  features.  Below  280  nm,  the  red  line  in  Fig.  7  corresponds  to  the  line-of- 
sight  spectrum  obtained  by  assuming  that  all  electronic  levels  of  NO,  except  the  A  and  B 
states,  are  in  LTE  at  the  measured  temperature.  The  NO  A  state  overpopulation,  pNo  a,  is 
determined  by  matching  the  height  of  the  NO  y  (0>1)  band  at  236.2  nm.  The  blue  curve 
(below  280  nm)  accounts  for  the  overpopulation  factors  of  the  NO  A,  C,  and  D  states, 
with  pNo  c  determined  from  the  NO  6  (0,5)  band  at  230.8  nm  and  pNo  d  assumed  to  be 
equal  to  pNoc-  The  NO  C  and  D  states  are  likely  to  be  in  equilibrium  with  each  other  due 
to  very  fast  collisional  exchange'^’^  owing  to  the  small  energy  difference  between  their 
potential  energies.  The  spectral  locations  of  relevant  NO  6  and  e  transitions  are  shown  in 
Fig.  8. 


200  205  210  215  220  225  230  235  240  245 

Wavelength  (nm) 


Figure  8.  Band  origins  for  NO  5  and  e.  The  spectral  emission  curve  is  the  non-LTE  NEQAIR2  calculation 
from  the  50  cm  test  section  (see  Fig.  7). 

Above  320  nm,  the  red  curve  corresponds  to  the  computed  equilibrium  spectrum, 
and  the  blue  curve  is  obtained  by  including  the  measured  pN2  c  in  the  simulation.  All 
measured  overpopulation  factors  are  summarized  in  Table  1. 


Table  1.  Measured  excited  state  overpopulation  factors 


Electronic 

state 

Transition 

name 

Term  energy 
(cm'‘) 

P 

40  cm 

P 

50  cm 

P 

65  cm 

NO  A 

NOy 

43966 

8 

500 

3.8x10"^ 

NOC 

NO  5 

52180 

25 

4600 

2.7x10^ 

NOD 

NO  8 

53085 

25 

4600 

2.7x10^ 

N2C 

2+ 

89137 

2.5x10^ 

7x10^ 

3.5x10" 

OH  A 

A^X 

32684 

1 

4.25 

18.3 
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Kinetic  Analysis 


The  air  plasma  and  nitrogen/argon  plasma  experiments  were  used  to  examine  the 
applicability  of  the  Dunn  and  Kang,i  Gupta  et  al. ,2  and  Park^-^  mechanisms  to  recombin¬ 
ing  atmospheric  pressure  plasmas.  As  it  was  concluded  that  the  Park  mechanism  most 
adequately  described  our  measurements  in  pure  air,9d0  only  this  mechanism  will  be  con¬ 
sidered  here.  The  kinetics  solver  CHEMKlNi  ’  was  used  for  one-dimensional  modeling 
of  the  plasma  chemistry  along  the  axis  of  the  flow.  Rates  to  describe  the  relevant  hydro¬ 
gen  and  argon  chemistry,  taken  from  Laux  et  al.,^  were  appended  to  the  Park  mechanism. 
Rates  for  and  N0'3.14  thermal  dissociation  by  argon  were  taken  from  the  NIST 
chemical  kinetics  database*^  (Table  2).  Dissociative  and  three-body  electron  attachment 
reactions  such  as: 

O2  +  e  <=>  O  +  O", 

O2  +  e  +  e  <=>  O2  +  e, 

while  certainly  important  in  low  temperature  air  plasmas,*^  are  not  considered  in  the  pre¬ 
sent  study  because  the  concentrations  of  the  negative  ions  are  negligible  at  temperatures 
above  2000  K.  Nor  did  we  consider  the  N4  dimer  formed  via  the  following  reaction: 

N2+  N+  +Mo  N^+M, 

as  it  is  likely  to  be  very  unstable  in  atmospheric  pressure  gases  at  temperatures  much 
above  room  temperature. 


Table  2.  Forward  rate  coefficients  in  Arrhenius  form:  A  T  Qxp{-E/T) 


Reaction 

A  (cmVmol/s) 

n 

^(K) 

O2  +  Ar  0  +  0  +  Ar 

1.2x10'^ 

0 

54280 

NO  +  Ar  N  +  0  +  Ar 

9.64x10^^ 

0 

74697 

To  calculate  centerline  velocities,  the  velocity  profiles  were  assumed  to  be  self¬ 
similar  to  the  measured  temperature  profiles  as  the  Prandtl  number  was  calculated  to  be 
near  unity.  The  velocity  profiles  were  then  scaled  to  match  the  4.7  g/s  mass  flow  rate  of 
argon/air/hydrogen.  The  inferred  centerline  axial  velocity  was  found  to  decrease  from 
1060  to  330  m/sec  between  0  and  65  cm.  The  resulting  temperature  history  is  displayed 
in  Fig.  9. 
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Figure  9.  Centerline  temperature  history  based  on  velocity  calculations  (air/argon  plasma). 

To  substantiate  the  velocity  predictions,  the  calculated  velocity  profiles  were  used 
with  the  measured  temperature  profiles  to  calculate  the  total  power  in  the  plasma  at  each 
axial  location.  The  calculated  decrease  in  plasma  power  as  a  function  of  distance  from 
the  nozzle  exit  was  then  compared  with  calorimetric  measurements  of  the  power  removed 
by  the  water-cooled  test  sections.  Satisfactory  agreement  is  found  (Fig.  10). 


Test  Section  Length  (cm) 

Figure  10.  Comparison  between  measured  and  calculated  power  removed  by  test  sections  in  the  air/argon 
plasma  experiments. 

The  predicted  nonequilibrium  factors  for  major  species  are  shown  in  Fig.  1 1 .  The 
mechanism  again  predicts  that  electrons  almost  exclusively  recombine  through  the  fast 
and  equilibrated  two  body  dissociative  recombination  NO^  +  e  N  +  O,  even  though 
NO"^  is  not  the  dominant  ion.  Therefore  the  three  body  recombination  reactions  control 
the  electron  overpopulation. 
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Figure  11.  Predicted  nonequilibrium  factors  for  the  air/argon  plasma. 


Despite  the  fact  that  overpopulation  factors  have  been  measured  for  various  mo¬ 
lecular  excited  states,  situations  exist  which  enable  the  data  to  be  used  to  evaluate  the  ac¬ 
curacy  of  the  chosen  ground  state  kinetic  mechanism. 


The  measured  NO  C  state  overpopulation,  for  instance,  can  be  compared  with  an 
upper  bound  value  corresponding  to  the  situation  where,  owing  to  strong  collisional  cou¬ 
pling,  the  C  state  is  in  Saha  equilibrium  with  the  free  electrons.  This  process  corresponds 
to  the  following  reaction  being  equilibrated: 

NO^  +  e  +  e  ^  NO  (C)  +  e. 

This  equilibrium  is  described  by  the  Saha  relation: 


NO'' 


^NO'^ 


^NOC 


Sno 


ItlnigkT 


3/2  ^vJ 


NO'' 


-E^IkT 


rvJ 

'NOC 


(3) 


where  Z''^  represents  the  rovibrational  partition  function  of  the  C  state  of  NO  or  the 
ground  state  of  NO’^,  and  is  the  energy  difference  between  the  NO  first  ionization 
potential  and  the  C  state  term  energy.  The  Saha  equation  can  be  expressed  succinctly  as: 

PnO  C  =  pNO+  Pe 

by  dividing  each  side  of  Eq.  3  by  equilibrium  values. 
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There  is  another  possible  mechanism  by  which  the  NO  C  state  population  can  be 
enhanced.  Indeed,  the  NO  C  state  predissociates  and  therefore  can  be  populated  through 
inverse  predissociation.i'^  Given  that  these  rates  are  fairly  large  when  compared  with  the 
rates  of  C  state  spontaneous  emission  and  collisional  excitation  processes, it  is  reason¬ 
able  to  assume  that  the  reaction 

NO  (C)  ^  N  (^S)  +  0(^P) 

is  approximately  equilibrated.  Since  this  reaction  involves  ground  state  N  and  O  atoms,  it 
follows  that: 

Pno  c  =  PnPo  • 

It  is  interesting  to  note  that  since  Park’s  mechanism  predicts  that  the  dissociative  recom¬ 
bination  reaction  NO"^  +  e  N  +  O  is  fast  and  equilibrated,  the  following  relation  holds 
at  all  locations: 

pNO+pe  2  PnPo 

Thus  both  Saha  coupling  and  predissociation  coupling  yield  the  same  predicted  upper 
limit  for  the  overpopulation  factor  of  the  NO  C  state.  Numerical  predictions  for  the  upper 
bound  on  pNoc  are  compared  with  measurements  in  Fig.  12a.  As  can  be  seen  in  Fig.  12a, 
the  upper  bound  determined  with  Park’s  mechanism  appears  consistent  with  the  measured 
overpopulation.  However,  since  it  is  expected  that  the  coupling  due  to  predissocia¬ 
tion/inverse  predissociation  is  very  strong  for  the  NO  C  state,  one  would  expect  a  closer 
agreement  between  the  measurements  and  the  predicted  upper  bound.  It  would  thus  ap¬ 
pear  that  Pn  is  too  high  as  po  is  relatively  close  to  unity.  Furthermore,  since  the  three 
body  recombination  of  NO  accounts  for  about  38  to  20  %  of  the  total  N  atom  consump¬ 
tion  between  15  and  60  cm  downstream  of  the  nozzle  exit  (the  balance  being  due  to  the 
Zeldovich  reactions  whose  rates  can  be  considered  well  known),  the  rate  of  N  +  O  +  M 
<->  NO  +  M  may  be  too  slow.  While  Park  concludes  that  this  rate  is  known  within  a  fac¬ 
tor  three, 3  examination  of  his  compilation  of  relevant  NO  thermal  dissociation  rate  meas¬ 
urements  seem  to  warrant  a  factor  ten  uncertainty  on  this  rate.  Thus,  Fig.  12b  includes 
predictions  based  on  calculations  in  which  the  third  body  efficiencies  for  this  reaction, 
excluding  M  =  Ar,  have  all  been  enhanced  by  a  factor  ten.  (The  argon  third  body  effi¬ 
ciency  was  not  enhanced  as  the  rate  of  thermal  dissociation  of  NO  by  argon  seems  to  be 
known  well. *5)  It  can  be  seen  from  Fig.  12  that  multiplying  the  third  body  efficiencies  by 
a  factor  10  improves  the  agreement  between  measurements  and  predictions  for  test- 
sections  <45  cm,  but  produces  a  recombination  rate  that  is  too  fast  at  longer  distances. 
Several  issues  must  be  examined  in  order  to  draw  more  definite  conclusions  about  these 
third  body  efficiencies.  First,  a  detailed  collisional-radiative  model  is  warranted  to  better 
determine  the  degree  of  coupling  between  the  NO  C  state  and  electrons  or  N  and  O  atoms 
via  the  two  processes  previously  mentioned.  And  second,  possible  departures  from  a 
Maxwellian  distribution  for  the  free-electrons  and  their  consequences  on  rate  coefficients 
should  be  considered.  A  preliminary  analysis  of  these  effects  (see  Appendix  B)  indicates 
that  departures  from  a  Maxwellian  distribution  may  occur  at  locations  >  35  cm.  These 
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various  issues  are  currently  being  examined  in  order  to  provide  a  better  assessment  of  the 
rate  of  NO  thermal  dissociation. 


Distance  from  Nozzle  Exit  (cm) 


Figure  12.  Comparison  between  calculated  and  measured  NO  C  state  nonequilibrium  factors.  Top  frame 
(12a)  uses  Park’s  mechanism.  Bottom  frame  (12b)  uses  Park’s  mechanism  with  a  rate  of  NO  thermal  disso¬ 
ciation  enhanced  by  a  factor  10. 

Similar  arguments  can  be  made  with  regards  to  the  interpretation  of  the  measured 
N2  C  state  overpopulation.  If  the  N2  C  state  is  in  Saha  equilibrium  with  the  free  electrons, 
we  then  have: 

PN2  C  =PN+Pe  • 

This  equation  defines  an  upper  bound  for  the  N2  C  state  overpopulation  corresponding  to 
full  collisional  coupling  with  free  electrons.  Another  upper  bound  corresponds  to  equilib¬ 
rium  between  the  N2  C  state  and  ground  state  nitrogen  atoms  as  a  result  of  predissociation 

and  its  reverse.  This  upper  bound  is  defined  by: 

_  2 

PN2C  -  Pn  • 

It  should  be  noted  that  this  is  the  model  employed  by  the  NEQAIR  collisional  radiative 
modeE^  for  the  N2  C  state.  Figure  13a  shows  the  comparison  between  the  measured  N2 
C  state  overpopulation  and  the  upper  bounds  just  described.  As  can  be  seen  in  Fig.  13a, 
these  two  upper  bounds  do  not  coincide.  This  is  because  Park’s  mechanism  predicts  that 
the  dissociative  recombination  of  Nj  is  not  equilibrated.  As  can  also  be  seen  in  Fig. 
13a,  the  measured  N2  C  state  overpopulation  is  larger  than  the  “upper  bound”  curve  cor¬ 
responding  to  predissociation  coupling.  This  may  indicate  that  Saha  coupling  is  stronger 
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than  predissociation  coupling  for  the  C  state  of  Na.  (It  should  be  noted  that  unlike  the  NO 
C  state,  not  all  rovibrational  levels  of  the  Na  C  state  predissociate  and  therefore  predisso¬ 
ciation  coupling  may  not  be  as  strong  as  in  the  case  of  NO  C.)  Nevertheless,  the  upper 
bound  defined  by  the  Saha  coupling  appears  consistent  with  the  measurements.  The  ef¬ 
fect  of  increasing  the  third  body  efficiencies  of  N  +  O  +  M  NO  +  M  by  10  is  shown  in 
Fig.  13b.  As  was  the  case  for  the  NO  C  state  analysis,  better  agreement  is  obtained  be¬ 
tween  the  Saha  coupling  curve  and  the  measured  Na  C  overpopulation  at  locations  <  45 
cm.  Yet  it  is  clear  that  a  more  detailed  collisional-radiative  model  is  required  to  interpret 
the  nonequilibrium  data. 


Figure  13.  Comparison  between  calculated  and  measured  Na  C  state  nonequilibrium  factors.  Top  frame 
(21a)  uses  Park’s  mechanism.  Bottom  frame  (21b)  uses  Park’s  mechanism  with  rate  of  NO  thermal  disso- 
ciation  enhanced  by  a  factor  10. 

Finally  we  compare  the  electron  overpopulation  factors  predicted  using  Park’s 
mechanism  with  the  measured  electron  overpopulation  factors.  At  0,  10,  and  15  cm,  di¬ 
rect  measurements  of  pe  were  made  via  the  Stark  broadening  of  the  Hp  line  at  486.1  nm. 
Error  bars  on  these  data  points  (Fig.  14)  reflect  uncertainties  on  measured  electron  num¬ 
ber  densities  and  on  centerline  temperatures.  At  the  exit  of  longer  test-sections,  the  inten¬ 
sity  of  the  Hp  line  was  very  weak,  and  therefore  no  direct  measurements  of  electron  num¬ 
ber  densities  could  be  made.  However,  Park’s  model  predicts  that  the  NO^  dissociative 
recombination  is  equilibrated  and  that  within  a  factor  3,  =Pe .  It  follows  that  pNO  c 
should  approximately  equal  (pe)^.  Figure  14a  compares  the  predicted  electron  overpopu¬ 
lation  factors  and  the  square  root  of  the  measured  Pno  c-  Figure  14b  shows  the  Pe  predic- 
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tions  obtained  if  the  third  body  efficiencies  mentioned  previously  are  multiplied  by  10. 
This  latter  curve  stops  after  42  cm  as  the  NO'^  +  e  N  +  O  reaction  is  no  longer  equili¬ 
brated.  These  results  support  the  same  conclusions  as  those  drawn  from  analyzing  the 
NO  C  and  N2  C  state  overpopulations. 


Distance  from  Nozzle  Exit  (cm) 

Figure  14.  Comparison  between  calculated  and  measured  electron  nonequilibrium  factors.  Top  frame 
(14a)  uses  Park’s  mechanism.  Bottom  frame  (14b)  uses  Park’s  mechanism  with  rate  of  NO  thermal  disso¬ 
ciation  enhanced  by  a  factor  10. 

Conclusion 

An  analysis  of  the  kinetics  of  recombining  air  plasmas  was  presented.  This  analy¬ 
sis  shows  that  electron  recombination  occurs  primarily  through  two-body  dissociative  re¬ 
combination  reactions.  As  a  result,  overall  electron  recombination  rates  are  significantly 
greater  than  those  for  three-body  eleetron  recombination,  and  the  extent  of  ionizational 
nonequilibrium  is  governed  by  the  rates  of  relatively  slow  three-body  atom  reeombination 
reactions  such  as  N  +  O  +  M  -»  NO  +  M  and  2N  -l-  M  -4  N2  -l-  M.  In  order  to  test  the 
predietions  of  three  reaction  mechanisms  widely  used  for  air  plasma  chemistry,  experi¬ 
ments  were  conducted  with  air,  nitrogen,  and  air/argon  plasmas  in  which  various  degrees 
of  chemical  and  ionizational  nonequilibrium  were  produced.  In  all  three  cases,  measure¬ 
ments  of  electron  number  densities  and  excited  electronic  state  concentrations  were  made 
in  order  to  assess  the  recombination  models  and  the  rates  of  the  controlling  three-body 
neutral  recombination  reactions.  Predictions  based  on  the  model  proposed  by  Park  were 
found  to  be  qualitatively  consistent  with  our  observations.  Before  more  quantitative  con¬ 
clusions  can  be  drawn,  several  issues  must  be  addressed.  First,  a  detailed  collisional- 
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radiative  model  must  be  developed  to  better  assess  the  degree  of  collisional  coupling 
between  electronic  excited  states  and  electrons  or  atoms.  Second,  the  simple  one¬ 
dimensional  model  used  to  predict  the  time-temperature  history  of  the  recombining 
plasma  should  be  refined.  Third,  the  consequences  on  the  rates  of  possible  departures 
from  a  Maxwellian  electron  energy  distribution  must  be  examined. 


Nevertheless,  it  already  appears  from  the  present  work  that  all  measurements  sup¬ 
port  the  proposed  air  plasma  recombination  mechanism.  In  addition,  while  the  rate  of 
three-body  N2  recombination  appears  to  be  accurate  within  a  factor  3,  our  measurements 
support  a  rate  of  three-body  NO  recombination  that  is  ~10  times  faster  than  the  rate  pro¬ 
posed  by  Park  and  Gupta  et  al.,  and  -100  times  faster  than  the  rate  of  Dunn  and  Kang. 


Appendix  A.  Argon  Continuum  Radiation 

The  continuum  radiation  resulting  from  electron-ion  free-bound  (e^),  electron- 
neutral  free-free  ( e“ ),  and  electron-ion  free-free  ( )  interactions  can  in  fact  be  esti¬ 
mated.  These  three  continua  correspond  to  the  following  processes: 

+  e  <-4  X  +  hv, 

X-i-e<-^X  +  e-hhv, 

X^  +  e  X^  +  e  -I-  hv, 

where  X  represents  Ar,  O,  N,  O2,  N2,  or  NO.  Argon  comprises  ninety  percent  of  the 
plasma  in  the  present  experiments  and  is  likely  to  be  the  major  source  of  the  observed 
continuum.  The  argon  continuum  emission  is  modeled  as  follows: 


£  =  £^*  +  £^”  +£^* 
ei  r'  n  ^  he 

e  ^  ig 

pi  he 

he 

e  ^ 

where  C\  =  l.bBlxlO''^^  Wm'^K'^^sr  ',  C2  =  1.026x10’^'',  ««,  and  n,-  are  the  electron,  Ar, 

and  Ar"^  densities,  respectively,  g,  is  the  electronic  ground  state  degeneracy  of  Ar"^,  the 
free-bound  Biberman  factor,20  Z,  the  Ar"^  partition  function,  the  free-free  Biberman 
factor,  and  Q(Te)  the  averaged  electron-neutral  collision  cross  section.  The  values  of 
^  and  Q(Te)  are  specified  by  Gordon.21  The  results  of  this  model  are  shown  in  Fig.  A.l 
for  a  2  mm  diameter,  uniform  temperature,  LTE  argon  plasma.  The  calculated  values  of 
2.5,  2.4,  and  2.2  |iW/cm^/sr  at  360  nm  are  consistent  with  the  constants  added  to  the  nu¬ 
merical  spectrum  within  a  factor  20.  Estimates  of  the  contribution  from  air  species  to  the 
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recombination  continuum  are  difficult  to  obtain  as  we  are  not  aware  of  calculations  and 
measurements  for  these  quantities  at  temperatures  below  9000  K.22.23 


250  300  350  400  450 


Wavelength  (nm) 

Figure  Al.  Total  radiative  recombination  emission  coefficient  for  a  1  atm  LTE  argon  plasma.  Calculations 
are  for  a  2  mm  diameter  plasma  observed  through  a  monochromator  with  a  0.44  nm  instrument  broadening. 


jendix  B.  Electron  Energy  Distribution 


Both  the  predicted  and  measured  chemical/ionizational  nonequilibrium  warrant 
consideration  of  the  validity  of  applying  Park’s  reaction  mechanism  to  this  data  analysis. 
In  particular,  it  is  possible  that  the  rates  of  importance  may  be  dramatically  altered  by  de¬ 
partures  from  equilibrium  in  the  Maxwellian  electron  energy  distribution  upon  which 
Park’s  rates  rest.  Sufficient  conditions  for  the  existence  of  a  Maxwellian  distribution  can 
be  obtained  by  consideration  of  the  first  Cartesian-tensor  equation  formed  from  the  elec¬ 
tron  Boltzmann  equation. 24  In  this  text,  the  authors  argue  that  the  electron-electron  colli¬ 
sion  integral  must  be  the  dominant  term  in  the  electron  Boltzmann  equation  if  a  Maxwel¬ 
lian  speed  distribution  is  to  exist.  An  order  of  magnitude  analysis  is  performed  and  the 
following  partial  list  of  conditions  is  derived: 


Vc/» 

h  ee 


«1 


(Bl) 


« 1  (B2) 

VggX 
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Cg  is  the  average  electron  speed;  Ue  the  electron  diffusion  velocity;  Vg/j  and  Vgg  the  elec¬ 
tron-heavy  and  electron-electron  collision  frequencies,  and  t  and  L  the  characteristic  time 
and  length  scales  for  macroscopic  change.  To  evaluate  these  expressions,  various  as¬ 
sumptions  must  be  made.  It  is  assumed  that  the  electron  energy  distribution  is  Maxwel¬ 
lian,  so  that  the  average  electron  speed  can  be  calculated  as  (8kT/7tme)‘'^".  T  is  assumed  to 
decrease  from  8000  to  2500  K  linearly  over  65  cm  such  that  Az/AT  =  1.2x10  '^  m/K, 
where  z  represents  the  distance  from  the  nozzle  exit.  A  100  K  temperature  drop  therefore 
occurs  over  1.18x10'^  m.  The  product  (ngT)  is  evaluated  at  LTE  conditions.  Finally 
V(neT),  evaluated  over  a  drop  in  temperature  of  100  K,  can  be  approximated  as  A(neT)IAz 
where  Az  =  1.2x10'^  m.  The  collision  frequencies  are  evaluated  with  equilibrium  number 
densities  and  the  electron  collision  cross  sections  found  in  Gupta.2 

The  first  of  these  conditions  results  from  the  requirement  that  the  electron- 
electron  collision  integral  be  much  larger  than  the  electron-heavy  collision  integral.  The 
next  two  expressions  assure  that  the  plasma  is  collision-dominated  with  respect  to  the 
electron  gas.  They  state  that  the  time  between  electron-electron  collisions  is  much  less 
than  the  characteristic  time  scale  for  macroscopic  change  and  that  the  mean  free  path  for 
electron-electron  collision  is  much  smaller  the  characteristic  length  scale  of  macroscopic 
change,  x  and  L  are  taken  as  the  time  to  traverse  the  longest  test  section  (1.29  ms)  and  the 
length  of  this  test  section  (0.65  m),  respectively.  The  final  condition  requires  that  the 
term  in  the  energy  equation  involving  the  electric  field  E’  is  small  compared  to  the  elec¬ 
tron-electron  collision  integral.  Numerical  evaluations  of  these  conditions  are  shown  in 
Fig.Bl. 
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Figure  Bl.  Evaluation  of  sufficient  conditions  for  the  existence  of  a  Maxwellian  electron  energy  distribu¬ 
tion. 

For  temperatures  less  than  -3800  K,  it  appears  that  these  sufficient  conditions 
may  be  violated  for  an  LTE  air/argon  plasma.  It  is  possible  that  the  observed  non-LTE 
behavior  of  the  air/argon  experiments  will  delay  a  departure  of  the  electron  energy  distri¬ 
bution  from  a  Maxwellian.  For  example,  the  inferred  pe  of  -500  at  65  cm  (Fig.  22)  would 
increase  Vge  and  thereby  definitely  prolong  the  validity  of  Eqs.  B2  and  B3,  and  possibly 
Eqs.  Bl  and  B4  depending  on  the  degree  of  chemical  nonequilibrium.  Considering  both 


20 


the  approximate  nature  of  these  calculations  and  the  fact  that  these  conditions  are  suffi¬ 
cient  but  not  necessary,  a  closer  examination  of  the  existence  or  lack  of  a  Maxwellian 
distribution  is  justified. 


3.3.  Significance  of  the  present  study 

The  results  from  this  work  have  important  implications  regarding  the  transport  of 
air  plasmas  and  the  development  of  methods  for  creating  and  sustaining  elevated  non¬ 
equilibrium  electron  number  densities  in  low  temperature,  atmospheric  pressure  air.  As 
discussed  in  this  report,  the  overall  recombination  of  electrons  results  from  and  is  con¬ 
trolled  by  the  recombination  of  neutral  species.  It  was  already  known  from  our  initial 
studies  using  existing  kinetic  mechanisms  that  the  three-body  recombination  of  NO  was 
much  faster,  by  about  a  factor  10,  than  the  recombination  of  N2.  From  the  present  work, 
it  appears  that  NO  recombination  may  be  even  faster,  by  about  a  factor  10,  than  previ¬ 
ously  thought.  It  is  clear  therefore  that  in  order  to  maintain  elevated  electron  number 
densities,  the  recombination  of  NO  must  be  inhibited.  In  future  work,  particularly  in  the 
framework  of  the  newly  started  Air  Plasma  Ramparts  program,  one  of  our  goals  will  be  to 
explore  ways  to  efficiently  sustain  air  plasmas  with  highly  dissociated  NO. 
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